The binary black-hole dynamics at the third-and-a-half post-Newtonian order in the 

ADM-formalism 
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We specialize the radiation-reaction part of the Arnowitt-Deser-Misner (ADM) Hamiltonian for 
many non-spinning point-like bodies, calculated by Jaranowski and Schafer to third-and-a-half 
post-Newtonian approximation to general relativity, to binary systems. This Hamiltonian is used for 
the computation of the instantaneous gravitational energy loss of a binary to 1PN reactive order. We 
also derive the equations of motion, which include PN reactive terms via Hamiltonian and Euler- 
Lagrangian approaches. The results are consistent with the expressions for reactive acceleration 
provided by Iyer- Will formalism in Ref. Q in a general class of gauges. 
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I. INTRODUCTION 

Compact binaries are the most promising sources of 
detectable gravitational waves for ground and space 
based interferometric experiments, TAMA300, GEO600, 
VIRGO, LIGO (Laser Interferometer Gravitational Wave 
Observatory) and LISA (Laser Interferometer Space An- 
tenna). Because of gravitational wave emission, the bi- 
nary separation decreases adiabatically in time while the 
frequency increases according to the relativistic version 
of the Kepler law. This process ends as soon as the first 
instable orbit is reached before the coalescence. Dur- 
ing the inspiraling phase the binary system produces a 
most likely detectable signal. The construction of wave 
patterns for data analysis requires the knowledge of the 
theoretical wave form to very high post-Newtonian or- 
ders [see, e.g., Rcfs. 3, 4]]. This assumes a good under- 
standing of the compact-binary dynamics. The secular 
effects governing the long-time evolution are of particu- 
lar importance. They result from the reaction force due 
to the quadrupole damping arising at the 2.5 order in 
power of v 2 /c 2 , v representing a typical orbital velocity 
and c standing for the speed of light in vacuum. For 
this reason, they are often referred to as the 2.5 post- 
Newtonian (PN) damping. The measurements resulting 
from the observation of the Hulse- Taylor binary pulsar 
PSR 1913+16 show excellent agreement with the values 
predicted by Einstein theory M. 

A while ago, Iyer and Will [6j derived the equations of 
motion for two non-spinning point-like objects including 
the first dissipative terms appearing at the 2.5 and 3.5PN 
level, i.e. the quadrupole damping and its 1PN correc- 
tion. They proceeded by postulating a balance equation 
between the instantaneous flux of energy and angular 
momentum in the far zone on the one hand, and the in- 
stantaneous loss of energy and angular momentum in the 
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system's near zone on the other hand. In a later paper 0], 
they detailed their derivation and proposed an alterna- 
tive method, based on the specialization of the Blanchct's 
radiation- reaction potentials to binary systems, which led 
to the same results Q • More recently, the dynamics of a 
two-body system at the 3.5PN order was derived by Pati 
and Will in harmonic coordinates [g using the integra- 
tion method for the Einstein field equations developed in 
Ref. 0. 

This article deals with the 3.5PN gravitational damp- 
ing of compact binaries regarded as point-masses. We 
base our investigation on the Af-body point-mass Hamil- 
tonian pj calculated by Jaranowski and Schafer within 
the Arnowitt-Deser-Misner (ADM) Hamiltonian formal- 
ism of general relativity [Tf]l | , which has proved to be very 
efficient in determining the approximate 3.5PN dynamics 
0, IH 13 G3 A remarkable feature of the ADM for- 
malism at this level is the absence of asymptotic match- 
ing between the near field and the far field given in a 
near-zone and a far-zone-defined coordinate system re- 
spectively. We can stay in one single global coordinate 
system in all our calculations, like in Ref. Q . Let us men- 
tion that the ADM formalism does not provide directly 
the balance equations between instantaneous losses and 
fluxes, even when assuming quasi-stationarity in the ra- 
diation emission. Although the ADM formalism gives us 
the conserved total energy, the balance between the lost 
energy of the matter system and the energy flux in the 
wave zone is proved up to 3.5PN order only [see, e.g., 
Ref. 0]. 

The main purposes of this work are (i) to specialize the 
A-body 3.5PN ADM Hamiltonian of Ref. Q to N = 2, 
(ii) to deduce from it the energy loss of a two-body sys- 
tem in general orbits, and (Hi) to calculate the equations 
of motions in the Hamiltonian form and in the Eulcr- 
Lagrangian form in order to determine the radiation- 
reaction force. The expression obtained for the energy 
loss is different from that derived in Jaranowski and 
Schafer 0] and in Iyer and Will @ , but the time-averaged 
expressions coincide in the case of quasi-elliptic orbital 



2 



motion. The plan of the paper is to reduce the original 
Hamiltonian for two point-masses. Next, we perform a 
partial time differentiation to get the energy loss. On the 
other side, we compute the radiation-reaction force from 
the Euler-Lagrangian equations of motion. The corre- 
sponding energy loss turns out to be equivalent to the 
preceding one after averaging. Our results are compared 
with those of Iyer and Will and the accessory gauge coef- 
ficients are determined explicitly. Hadamard regulariza- 
tion is systematically employed. 

We use units in which 16irG = 1 = c, where G is 
the Newtonian gravitational constant, but G and c are 
restored in the presentation of our final results. In our 
notation, t is the coordinate time. Characters in bold, 
like x = (x l ) , represents a point in the 3-dimensional 
Euclidean space M 3 endowed with a standard Euclidean 
metric. The scalar product is denoted by a dot. Let- 
ters a, b, . . . are particle labels, so that x a S M 3 de- 
notes the position of the ath particle. We also define 
r a := x - x a , r a := |r Q |, n a := r a /r a ; and for a ^ b 
r 6 := x a - x 6 = r h - r a , r ab := \r ab \, n ab := r ab /r ab 
| ... | stands here for the Euclidean length of a vector. 
The momentum vector of the particle a with mass pa- 
rameter m a is denoted by p a = (p a ;)- A dot over a 
symbol, like in x a , represents its total time derivative. 
The partial differentiation with respect to x 1 is denoted 
by di, or equivalently by a comma, i.e. di<p = 4>.i- 



takes the form 

H = H [x a ,p Q ,/iT. T ,7r^ TT ] . (2.2) 

The Hamiltonian equations of motion for body a are 

9H dH , nn . 

Pa = -Tr— and x a = - — , (2.3) 



<9x a 



dp a 



and the field equations for the independent degrees of 
freedom read 



d 

5? 



r tjTT 



rTTij 



at ij ij 5ir klTT ' 



(2.4) 

where (6 . . .)/(S . . .) denotes the Frechet derivative [see, 
e.g., Ref. and where STT'^ is the TT-projection op- 
erator [see, e.g., Eq. (2.17) of [13]. 

The Hamiltonian H can be uniquely decomposed into 
a matter part, H mat , depending on the matter variables 
only, a field part, i? field , depending on the field variables 
only, and an interaction part, H mt , depending on both 
sets of variables, matter and field, in a way that it van- 
ishes if one set is put to zero. In short, the full content 
of the field-plus-matter dynamics at the 3.5PN order is 
included in the Hamiltonian 
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II. 3.5PN FIELD EQUATIONS 

The ADM formulation is based on a (3 + 1) splitting 
of space-time. The canonical variables are, on one hand, 
the projection of the 3-metric on the space-like hy- 
persurfaces x° — t = const, in a suitable coordinate grid 
and, on the other hand, its co njugate momentum ~k %3 . In 
the so-called ADM gauge 0, lL2t Tl5| , the spatial metric 

gij decomposes into a diagonal part (l + Sij, plus 
a part hj^ transverse and traceless with respect to the 
Euclidean metric Sij , so that 
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<f> + h 



TT 



(2.1) 



with djhj? — 0. In addition, the field momentum tt 4 - 7 
is traceless: 7r" = 0. Let us consider now hJ n T as the 
new field variable, and denote its conjugate momentum 
by 7r yTT . The degrees of freedom associated with -P- 7 = 
n v _ jrij'TT anc i the potential <j> are removed by solving 
the constraint equations. 

After inserting the values of ^ and into the gen- 
eral relativistic Hamiltonian of the isolated system under 
consideration in an asymptotically flat space-time, the 
latter simplifies to a surface integral often referred to as 
the reduced Hamiltonian H. It is a functional of the 
independent degrees of freedom hJ 3 T and tt' 1: > tt 01 ■ 

For a system of N point-like bodies with position vec- 
tors x a and momenta p (a = 1, . . . , N), the Hamiltonian 



The notation < 3.5PN indicates that all orders lower or 
equal than 3.5PN are included. In what follows, only 
the interaction part of the Hamiltonian l|2.5 r will be 
needed . It has been explicitly computed at the dissi- 
pative 3.5PN order in Ref. [l| by Jaranowski and Schafer. 



Its truncation at the 2.5PN level. H, 



N , takes the well- 
(2.6) 



known form presented in Refs. |12lll 

#2 n 5PN ( x a, Pa, t) = 5?T 0) X{i)ij ( x a, Pa) 

with 
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X(i)ij ( x a,Pa) 



1 
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E 



1 \ - \ - m a m b 
16^^^ 



Tab 



(2.7) 



The function X(4)ij (t) comes from the post- Newtonian ex- 
pansion of the retarded integral giving hJ 3 T as a function 
of x a and p a Q. As dH/dx a and dH/dp a are differ- 
entiated for constant hj? and 7r yTT , the positions and 
momenta appearing in Eq. (|2.7|l are not affected by the 
operators d/dx a and d/dp a - In order not to confuse the 
latter particle variables with the ones that are already 
present in the original Hamiltonian (|2.2|l and are affected 
by d/dx a and d/dp a , we shall mark them with a prime 
symbol: x tt /, p a > . Thus, in our notation, XU)ij{t) denotes 
the time derivative of X(4)ij( x a': Pa')- 
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The 3.5PN Hamiltonian fl| n | PN reads 1 



H 3 n 5PN ( x a,Pa,*) = 5n X(A)ij (*a,Pa) 



n Uj (t) + u 2ij (t) + n 3 y (0 



57T 



X(4)ij (*) Illy (X Q , Pa) + n 2 y (X a , 



57T X(4)ij (*) n 3-y (x a , Pa) + X(4)ij (*) Q'ij (x a , Pa, t) + Q" tJ (x a , t) 
Q3 



dt 3 



R' (x a ,p a ,t) + R" (x Q ,t) 



(2.8) 



_ffg n g PN is parametrized by the following functions 
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(2.9b) 
(2.9c) 



1 We observe that there are few misprints in the expressions for 
ffg n gp N , appearing in [lj. We make appropriate changes in those 

expressions and their positions are marked by | . . . | . 
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(2.9d) 



with s aa 'b' := r aa / +r ao / +r /(,/ . Variables originating from /i? T and 7r* : ' TT are primed. It is indeed crucial to distinguish 
the particle positions and momenta inside and outside the transverse-traceless (TT) quantities at this stage. 



The other four functions which enter H™^ PN are 
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with s a6a > := r ab + r aa > + r ba , 
2 / 1 

i? (x a ,p a ,i) := 
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Now, we are in a position to reduce the ./V-body 3.5PN 
Hamiltonian (|2~%|) to N = 2. The result will serve as a 
base for the derivation of the gravitational energy loss of a 
compact binary moving in general orbits, like in Ref. , 
and provide a reference to check the equations of mo- 
tion that will be calculated afterward from this two-body 
Hamiltonian. 



relates to their behaviors with respect to time reversal: 
#<3P N is symmetric whereas -ff<3 S 5PN is antisymmetric. 

The total time derivative of i?o.5PN is equal to its par- 
tial time derivative, the other terms canceling by virtue 
of the equations of motion. Only the dissipative part of 
-ff<3.5PN depends explicitly on time, through the substi- 
tution of the transverse-traceless variables that have been 
replaced by their actual values. We get therefore 



III. ENERGY LOSS OF A TWO-BODY SYSTEM 

Let us denote by i?o.5PN the Hamiltonian resulting 
from -ff<3.5PN, Eq. l|2.5[) after writing the field part as a 
function of x a , p a , x a , p Q and subtracting its 2.5PN and 



3.5PN contribution H$ e $ N and H$%% 



By definition, we have 



from i/| c 3 ld 5PN . 



H 



<3.5PN 



Rinat i rrficld 

-"^.SPN ' -"O.SPN 

i rrint 
"+" -"O.SPN ) 



H. 



field 
2.5PN 



H. 



field 
3.5PN 



(3.1) 



which can be split into a conservative and a dissipative 
part 

-ff<3.5PN (x a , Pa, X a , p a , t) =7J<3 1 p N (x a , p a , X a , p a ) 

+ H<3%PN ( X a, Pa, t) , 

(3.2) 



where 
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(3.3a) 
(3.3b) 



The conservative Hamilton function corresponding to the 
expression (|3.3a|) is given in Refs. |18lll9| . With the help 
of the coordinate transformation described in Ref. |13j . 
it is possible to eliminate the dependence on x a and p a . 
The difference between H%>™ n and -ff<3 SS 5PN in Eq. (E3 

2 All algebraic computations arc performed using the commercial 
software Mathematica 20) • 



dt 



-ff<3.5PN (x a , p a , X a , p a , t) 

d 



TT7-ff<3.5PN (x a , Pa,X a , pa, t) 
Ot ~ 



9t s 



<3.5PN 



(x a ,p Q ,i) 



(3.4) 



The instantaneous energy loss of the matter system due 
to the gravitational wave emission is defined as 



£<3.5PN (*) '■= -^--H<3.5PN (x a , Pa, X a , pa, t) . (3.5) 

Out task is now to compute £<3.5pn for a system of 
two point-like bodies with masses mi and m 2 . 2 We start 
by rewriting the right-hand side of Eq. I|3.5|l [by means 
of Eqs. (|3.3bjl and (|3~^|) ] in the form 
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d_ 

dt 
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rrint \ 



(3.6) 



Next, we insert the formulas 1|2.6[) and (|2.8I) into Eq. (|3.6[l . 
and specialize all sums of Eqs. (|2.7|) and l|2.9|l for N = 2. 

Notice that the triple sum X) S S disappears in rela- 



a b^a c^a,b 



tions H2.9afl and 12.9b|) . For clarity, we also reintroduce 
the Newtonian gravitational constant G and the speed of 
light c. 

In differentiating the Hamiltonians l|2.6|l and (|2.8|l with 
respect to time, we replace the time derivatives of the 
primed coordinates and momenta according to the 1PN 
equations of motion: 
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(3.7a) 
(3.7b) 
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Gniy 7712 



m 2 ' 2 o TO l' 2,7/ \ . q / \ / 

3 p v - 3 p 2 , + 7 (pi/ • p2'J + 3 (ni/ 2 ' ■ Pi') (ni>2> • p 2 ' 



m 2 ' 

"'2 



n V2 > 
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(3.7d) 



The calculation of the third time derivative of R' and R" is particularly tedious. 

The particle momenta p a / a < is expressed in terms of the particle coordinate velocities v a /a' U P to the 1PN order 
[see, e.g., Eq. (4.1) of H3] by: 



Gmi m 2 

2c^ 2 ri2 

1 , Gm\rri2 
p 2 = m 2 v 2 + 7r^-m 2 v 2 V2 + — -= 

2c z 2 c z ri2 



6vi - 7v 2 - (ni 2 ■ v 2 ) ni 2 
6v 2 - 7vi - (ni2 • vi) ni2 



(3.8a) 
(3.8b) 



Analogous relations hold for pi' and p 2 '. 

The time derivatives of the primed coordinates Xi' and x 2 ' coming from the differentiation of -ffo.spN in Eq. I|3.5(l 
equate the primed coordinate velocities Vi' and V2'. We eliminate all accelerations by making use of the 1PN equations 
of motion [see, e.g., Eq. (1.5) of fl5|]: 
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(3.9b) 



At this stage, we can identify positions and velocities of 
particles inside and outside the transverse-traceless vari- 
ables [i.e. the primed and unprimed quantities]. As the 
limit Xi' — > Xi, x 2 ' —* X2 is singular, we shall resort to 
the Hadamard regularization procedure described in Ap- 
pendix A [for a more detailed investigation, see paper |22j 
and the references therein]. The computations are ex- 
tremely long, but straightforward. 

To write the final formula for the energy loss in a more 
compact way, we introduce the total mass of the system 
M := mi +i7i2, the reduced mass /i :— mi m 2 /M, as well 
as the ratio parameter v :— [ijM. The individual masses 
mi and m 2 are given in terms of /i and v by: 



m 2 



(3.10b) 



assuming mi > m 2 . We also express the velocities of the 
bodies in the center-of-mass frame as a function of their 
relative velocity v = vi 2 := vi — v 2 and of their sepa- 
ration r = ri2 := xi — X2 — Y2 — i"i with v = f . This is 
achieved by means of the relations linking the positions 
of the individual particles to r at the 1PN order [see, e.g., 
Eq. (2.4) of H3]: 3 



xi 



mi 



fi{mi - m 2 ) 
2c 2 M 2 



GM 



(3.11a) 



mi 



(3.10a) 



3 At the 1PN order the Eqs. l3Jlt in ADM gauge and in harmonic 
gauge are identical. 



7 



x 2 = 



(i (J-(rni — m-i) I 2 GM 
^ + 2c 2 M 2 ' V 



(3.11b) 



where r = |r|. By differentiating Eqs. (|3.11() with respect 
to time and eliminating the temporal derivatives with the 
aid of the Newtonian equation of motion 



GM 



-n, 



(3.12) 



where n := r/r, we are led to: 
2vv VI - 4^ 



vi = 



l + x/l - 4^ 



2c 2 



(v + (n • v) n 



(3.13a) 



-2^v VI - 4z/ 

V2= l-Vl-4, + ^^ X 
2 __ G/i 



v ( v + (n • v) n 



(3.13b) 



Finally, we arrive at the instantaneous energy loss [24 



(i) 

<3.5PN 



4 G 2 M 3 z/ 2 LG 2 M 2 



15 c 5 r 



1 G 2 M 3 v 2 
105 c 7 r 3 



+ 



(5 - 104^) 



llv 2 -9(n-v) 

G 3 M 3 



GM 



llv 4 - 60 (n • v) 2 v 2 + 45 (n • v)" 

G 2 M 2 



(-1615 + 776i/)v 2 + (2733 + 92z/) (n • v) z 



- (1308 + 1372j/)v 4 + 4(4158 + 1577z/) (n • v) 2 v 2 - 3(5174 + U96v) (n • v) 



GM 



r 



(164 + 1064^)v B + 3(716 + 1562z/) (n ■ v) 2 v 4 - 15(324 + 56^) (n ■ v) 4 v 



+ 105(28- 22;/) (n ■ v) 



(3.14) 



The comparison of the expression above for the instan- 
taneous gravitational energy loss with the known instan- 
taneous far-zone flux [see Eq. (51) of 25] corrected by an 
erratum and equation below Eq. (3.40) of 26]] shows that 
they are not either identical. This is nothing to worry 
about since the form of the instantaneous gravitational 
energy loss in an isolated system is as well coordinate as 
representation dependent. 

We also observe that the expression (|3.14|l is differ- 
ent from formula (73) [in the work of Jaranowski and 
Schafer jlj]] at the 3.5PN order. However, the difference 
may be expressed as a total time derivative. To verify 
that, we computed (£<3 5PN ) and (£< s 3 5PN ), the orbital 



average of 5PN and 



£<3.rpn respectively, using the 
prescription given in Ref. |26| . for quasi-elliptical orbital 
motion and found (£< s 3 5PN ) = (£<3 5pn)- So, we recover 
Eqs. (4.20) and (4.21) of Ref. |2fj7 The important point 
is that all averaged losses coincide. 

In the following, we shall make the assumption that 
the binary has entered the inspiraling phase, so that its 
relative motion is quasi-circular. The investigation ap- 
pears to be simpler in this case, since all terms involving 
the factor (n • v) vanish. The averaging is performed 



by expressing C 



i) 

<3.5PN 



as a function of E. According 



to equation (3.35) in Ref. |26J, we have for the energy 



(reduced by the factor n) 



GM 



2 r 

GM " 



(l-3z,)v 4 

8c 2 



2rc 2 



(3 + v)v z + v (n • v) 



GM 



(3.15) 



The right-hand side depends on the relative velocity v 
and the relative distance r, whereas the symbol E in the 
left-hand side may be viewed as an adiabatic parameter. 
By solving Eq. I|3.15|l for v 2 , we find 



= 2E 



2GM 



1 



£ 2 (-3 + 9i/) 



GM 



-£(-12 + 16v) 



G 2 M 2 



(-10 + 7*/) 



(3.16) 



Next, we substitute Eq. (|3.16|) to v 2 in the instantaneous 
energy loss (|3.14|l . and eliminate the relative distance r 
in favor of the reduced energy with the aid of Eq. (43) in 
Ref. [23: 



GM 
~2~eT 



1, ,E 



(3.17) 



8 



As a consequence, the time-averaged energy loss special- 
ized to the case of circular orbits (c.o.) reads 



<3.5PN 



1024 v 2 {-Ef 
~ c 5 G 



(-E) 



13 

168 



(3.18) 



We recover formula (4.20) in Ref. [2(j with eR = (due 
to the fact that the orbits are quasi-circular). This is in 
agreement with Eq. (81) of Ref. [25j . Processing Eq. (|4.6|) 
or Eq. (73) in Ref. |l| in a similar way ends in the same 
result. 

The instantaneous gravitational energy loss £<3.spn 
given in Eq. I|3.14|l will be taken as reference to check the 
equations of motion that will be determined in the next 
section. 



IV. HAMILTONIAN EQUATIONS OF MOTION 

In this section we compute the reaction part of the 
equations of motion 



Pa 



url <3.5PN 



and 



C,J2 <3.5PN 



dp a 



(4.1) 



in a fully explicit form in terms of xi,X2,pi and p2. 
For this goal, we have to differentiate the Hamiltoni- 
ans l|2.6|) and (|2.8|l with respect to the particle coordi- 
nates x a and the particle momenta p a with a = 1, 2 and 
apply the Hadamard regularization procedure described 
in Appendix A. We find j2J] : 



1 G 2 
x i =-z< -o 



24 / w x TO 2 /24 . , 2 16 2 \ 16 , 

— (n ■ pi) (n • p a ) + — I — (n- Pl ) - y Pl 1 + — ( Pl ■ p 2 ) 



8m 2 
3m! 



(n- P i) 



-15 (n ' P2) 



16 . 

Pi + y (n • Pi)p 2 



1 G 2 



m,2 ( 36 
m\ \ 7 



\4 60 2 2 , 4 4 

(n-Pi) - — (n-pi) Pi + — Pl 



1 / 148. , 3 , , 748 , w , , 12 , , 2 , x 172 , , 
^2 — f ( n ' Pi) ( n ' Pa) + — (n • Pi) (n ■ p 2 ) p{ + — (n • Pl ) (p x ■ p 2 ) + — p{ ( Pi • l 

1 /104 2 , 2 324 2 2 432 136. 
-z-(n-pi) (n-p 2 ) - — (n-p 2 ) p 1 - — (n • Pi ) (n • p 2 ) ( Pi • p 2 ) - — (pi • 

m,\m,2 \ 7 35 35 35 



Pa)' 



132 , , 2 2 172 o 2 \ 1/8. w , 3 292 , . 2 , 
— (n • Pl ) p\ - -p\p\ j + (- (n • Pl ) (n . p 2 ) 3 + - (n • p 2 ) ( P i • p 2 ; 



m 2 / 324 3 1672 2 

— 3 -^r(n-Pi +-77^-(n-Pi)Pi 
mv V 35 105 



— (n • pi) (n • p 2 ) p\ + — (pi • p 2 ) p 2 

1 /64 , , 2/ , 1136 , ,0 184, 
^y(n-Px) (n • p 2 ) - — (n • p 2 ) Pl - — (n • Pl ) ( Pl • 

76 , , . ,2 376 , w , 992 , , 2 \ 
— (n ■ pi) (n • p 2 ) + — (n • p 2 ) (pi • p 2 ) - — (n ■ Pl ) p 2 J 



1 

'2 



l I 208 , ,3 112 , . 2 ' 

{"to (n ' P2) +l^( n 'P 2 )P 2/ 









): 


pi + 











1 / YD . 3 108 2 

-2 y(n-Pi) - — (n. Pl)Pl 



1 



1 /292 



992 , , 2 , x 664 , , , 400 , 
"35- ( n ' Pi) ( n ' P2) + yg- (n ■ p 2 ) p x + — (n ■ pi) (pi • p 2 ) 



1 /292, w x2 1928 , w 52 , . Al 1 

yi^^-Pi)^'^) -y^(n-p 2 )(Pi-p 2 ) + I ^(n. Pl )p 2 Jjp 2 } 

G 3 f \ml ( 7862 , s2 4129 A / 248 , , 2 3756 , 



91 

852 



m| 



7862 
105 



2 21 06 , 
35" Pl " IS" (Pl ' P2 



t-mi 



105 

/5032 
V 105 



(n • Pi) + yyPi + ™a ( — y (n • Pi) 



3756 
~35~ 



(n • pi) (n • p 2 ) 



, 158 , , 2 484 , x 97 
(n • Pl ) (n • p 2 ) - — (n • p 2 ) - — (pi • p 2 ) + — 



'617m 2 . 



i 2 . /2588. . 431 , A 64 

i5m7 (n ' Pl) + m2 Vm (n ' Pl) " (n ' P2) J - i05 mi (n ' P2) 



158 
~5 

64 
105 



Pi 



pi 



n 



9 



+ 



1469 , . / 484 , . 103 , 

-^-m 2 (n • pi) + mi I -— (n ■ pi) + — (n • p 2 ) 



P2 



G / 2 2 

— j- ym 1 m 2 + mim^ 



x 2 = (1 ^ 2) . 

For the time derivative of the particle momenta, we get 
1 fG 3 



—mim 2 (n • pi) - — m 1 m 2 (n • p 2 ) 



16 2 16 2 
n —m 1 m 2 pi + — m 1 m 2 p 2 



(4.2a) 
(4.2b) 



1 G 2 



m 2 
m? 



116 



10 (n • Pl )° + ^ (n • Pl ) 3 P? - ^ (n ■ Pi) pf 



198 



7 



35 



+" 32 ( 30 ( n ' Pi) 4 ( n ■ P2) - ~y ( n • Pi) 2 ( n ' Pa) Pi + y (n • P2) Pi - ^ (n • pi) 3 (pi • p 2 ) 



568 , , o , . 
+ -gg-(n-Pi)Pi (Pi -p 2 ) 



1 



mim 2 



30 (n • pi) J (n • p 2 ) z + 6 (n • pi) (n • p 2 ) 2 p 2 



+ 44 (n • pi) 2 (n • p 2 ) (pi ■ p 2 ) - ^ (n • p 2 ) p 2 (pi ■ p 2 ) - ^ (n • pi) (pi • p 2 ) 2 

35 35 



30 



+ ^(n-pi) 3 p 2 -^(n.pi)p 2 p 2 ) +£2 I 10(n.pi) z (n.p 2 ) J + ^(n.p 2 r P 2 



38 



7 



35 



1 



22 



\3„2 



7 



124 2 / \ . 208 ,2 30 2 / x 2 
— (n-pi)(n-p 2 ) (pi • p 2 ) + — (n- p 2 ) (pi • p 2 ) -y(n-pi) (n-p 2 )p 2 



124 

-2 (n • p 2 ) p 2 p 2 . + — (n • pi) (pi • p 2 ) p 2 



n 



m 2 ( 26 , 4 2 2 

— (-— n-pi 4 + 4 n-pi 2 p? 
m? V 7 



62 
105 



Pi 



1 / 40 . .3 . .16. .. . 8 . 2 1 \ 32 o / \ 
+ mf V~ ' ' ~ "5 ' P ' P Pl ~ 5 ' ' + 105 Pl ' 



H — ( (n • Pi) 2 (n • p 2 ) 2 + \ (n • p 2 ) 2 p 2 - | (n • pi) (n • p 2 ) (pi • p 2 ) 

mim 2 V 7 7 5 



16, ,2 2 , 2 2 22 2 2 \ , 1 f 12, w ,3 , 36, , 2 

In, • rin I -I I n • ni ' •-- 

21 



105' 



— (Pi ' P2) + 7; (n • pi) p 2 - — Pi P2 ) + ^2 \y — Y ( n ' Pi) ( n ' P2) + ^ ( n • P2) (Pi • P2) 



35 



4 4 x 

+- (n • pi) (n • p 2 ) P2 - — (pi • p 2 ) p| 









): 


pi + 








mf \ 



12 



22 



— j(--(n-pi) + — (n-pi) Pi-— Pi 



1 f 12 , <i, , 12 / w \ 2 204 ^2, , 28 2 . 

\ I — (n • Pi) (n ■ p 2 ) - — (n • pi) (n • p 2 ) p x - — (n • p x ) (p x • p 2 ) + — Pi (pi • p 2 ) 



mim 2 \ 7 
1 / (i 



35 



66 

+ Z?2 ( -f ( n ' PiT ( n ' P2) 2 - jjg (n • p 2 )^ Pi + ^ (n • Pl ) (n • p 2 ) ( Pl ■ p 2 ) - ^ ( Pl • p 2 ) z 



^2 + ±92 
35 



15 J 

40 



21 



6 , \2 2 ^22 

+^(n-pi) P 2 + gPiP 2 



G 3 



ml ( 52 3 
— =- (n- pi) 

mi V 5 



1408 
T05~ 



(n-pi)p 2 



20 



84 



124 



+ m 2 — (n-pi) + — (n-pi) (n • p 2 ) - — (n • pi) p : 



138 



7 



105 

2 



105 



(n- p 2 )p 2 



3128 . , . A / 20 , , 2 , . 464 , . 2 68, , 2 

— jqq- ( n ' Pi) (Pi ' P2) I + m 1 I — — (n • pi) (n-pa)-— (n-pi)(n-p 2 ) + — (n-p 2 ) Pl 

24 , , , . 2624 , , , . 1472 , , , 

+ (n • Pi) (Pi • p 2 ) + -j^r (n • p 2 ) (pi • p 2 ) + (n • pi) p^ 



m 2 /48 3 324 2 

+ m7(y (n ' P2) -^ (n ' P2)P2 



n 



m 2 / 62 2 82 2 



^ --(n-pir + — p 2 
mi V 21 105 



10 



m 2 



128 a 
T05 (n ' Pl) + T05 



1664 . w 8 2 188 . 

(n • pi) (n • p 2 ) + —Pi + — (pi • p 2 , 



-i" i I gg (n • Pi) (n • p 2 ) - ^ (n • p 2 ) 2 



8 / x 34 2 
log CPi ■ Pa) - iQgPa 



Pi 



+ 



134 j 
<■« < I "2j- ( n - Pi) 



254 
105 



Pi 



'148 , , 2 2024 . 

"" 1 Io5 (n ' Pl) -ioy( n -Pi)( n -p 2 ) 



52 
105 

174 



Pi 



292 
105 



(Pi ' P2 



'30. , 2 14 2 
— (n-pa) - — p 2 
m 2 V 7 5 



P2 



3574 
T05 



■m x m\ (n • p x ) 



5 

692 



raj ra 2 (n • p 2 J - 
, 3268 



3838 , , 3782 , 

( n -Pi) + -T7Tr ( n 'P2, 



105 



21 

P2=(1^2). 



105 



TOiTO 2 



Pi 



1116 



105 
164 



35 



-m\m2 



P2 



(4.3a) 
(4.3b) 



To check the correctness of the expressions for x a , p a , we employ them to compute the instantaneous energy loss 
of a binary, at 2.5PN and 3.5PN orders by differentiating 1PN accurate conservative Hamiltonian expression [see, 
Hamiltonian (2.5) of Ref. truncated at the 1PN level]. 4 The contribution to the Hamiltonian above comes purely 
from the matter part and it reads 



^<ipn( x oi Pa) 



p 2 p 2 , Gm\m,2 1 
2mi 2rri2 r 8c 2 

| ? (Pi - Pa) (n • Pi)(n ■ p 2 ) 



Pi 



p! 



Gm\rri2 


[-3 (4 - 


_P| 


2rc 2 


. \ m i 


mi 



TO1TO2 

Differentiation with respect to time leads to 

d 



tt mat 



pi rrmat 
a - rz <lPN . 



-xi 



TO1TO2 



pi rrmat 
U£1 <1PN ■ 



G 2 m 1 m 2 ( TO i + m 2 ) 
2r 2 c 2 



9x 2 



-x 2 



pi rrmat 
ua <1PN . 



9pi 



-Pi 



rrmat 
c ' rJ <lPN . 



C?P2 



-P2 



(4.4) 



(4.5) 



We then replace the time derivatives of x 1; x 2 , Pi and P2 occurring in Eq. (|4.5I) by their expressions l|4.2|) and (|4.3|1 . 
At last, we write the particle momenta in terms of the particle coordinate velocities with the help of relations (|3.8|l . 
We introduce the relative velocity, the relative position, the total mass and the parameter v, and multiply by — 1 
according to definition (|3.5|l . In this way, we obtain for the luminosity |24| : 



C 



(2) 

<3.5PN 



= ~dt H ^ 



<1PN 



4_GW 

15 c 5 r 3 

1 G 2 M 3 u 2 
+ 105 c 7 r 3 



24v 2 - 22 (n • v) 2 



GM 



105 



G 3 M 3 



(4709 + 384i/)v 2 + (4821 + 280i^) (n • v) 



G M 2 



(914 - 784i/)v 4 - 2(1313 - 1894i/) (n • v) 2 v 2 + 3(508 - 1012i/) (n • v) 4 
(62 - 344^)v 6 + 3(58 + 284^) (n ■ v) 2 v 4 - 450(3 - Av) (n ■ v) 4 v 2 

1050(l-2i/) (n-v) e 



GM 



(4.6) 



Though this result differs from the instantaneous energy 
loss (|3.14() . one will see by using the method presented 

4 The 2PN terms do not play any role in the computation of the 
2.5PN and 3.5PN energy loss, since they first couple with the 
reaction force at the (2 + 2.5)PN = 4.5PN order. 



in Ref. |26| for averaging that it leads to the same aver- 
aged value, in agreement with Eqs. (4.20) and (4.21) of 



11 



Rcf. 

V. EULER-LAGRANGIAN EQUATIONS OF 
MOTIONS 

The aim of the present section is to compute the PN 
reactive corrections to x Q employing two slightly differ- 
ent methods. We will see that the computation of these 
Euler-Lagrangian equations of motion does not involve a 
Lagrangian in cither of these methods. However, we will 
require PN accurate relations between particle momenta 
and its coordinate velocities. We turn to derive these 
relations in the next subsection. 



A. Legendre transformation 

For the present computation, wc require 1PN, 2.5PN 
and 3.5PN corrections to the familiar p a = rua^a rela- 
tions. So, we have to invert the equations (14.21) giving 
xi = /(xi, X2, pi, P2) and x 2 respectively, in order to 
determine pi = /(xi,x 2 ,xi,x 2 ) = /(xi,x 2 , Vi, v 2 ) and 
P2. This will provide the link between the particle mo- 



menta and the particle coordinate velocities up to the 
order 1/c 7 . 

As the first step, we add Newtonian and 1PN contri- 
butions to x a to the right-hand side of Eqs. 1)4. 2|) . These 
corrections are obtained from Eqs. (|3.7a[l and (|3.7b|) after 
replacing primed quantities by their unprimed counter- 
parts. Through the expression presented below refer to 
the first body, it should be noted that similar expressions 
hold true for the second body. It is the full coupled sys- 
tem that must be inverted. The time derivative of the 
position Xi has the form 

xi = vi = — + \a 2 + -^rA 5 + —A 7 , (5.1) 

mi cr c° c' 

whereas the structure of pi after inversion is 

111 , . 

pi = mivi + -^B 2 + —B 5 + —B 7 . 5.2) 

c z c° c' 

The coefficients B 2 , -B5 and Bj in Eq. (|5.2() are computed 
iteratively in terms of A 2l A§ and A-j of by means of 
Eq. 1)5.1(1 . The first-order coefficient B 2 is read directly 
from Eq. (|3~%a|l . We get 0]: 
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(5.3a) 



(5.3b) 



B. First method to calculate the Euler-Lagrangian equations of motions 

The calculation of the accelerations makes the connection with the Euler-Lagrangian approach. It is achieved in a 
first stage by evaluating the total time derivatives of the functions xi = /(xi, x 2 , pi, p 2 ) and x 2 = /(xi, x 2 , pi, p 2 ) as 
given by the Eqs. (|4.2I) . It is worth stressing that this method does not need the explicit knowledge of the Lagrangian. 
We can write indeed: 



d_ 

dt 
d 



<9xi . 9xi . 9xi . 9xi 

Xi =3T X 1 = n — x l + a — x 2 + a — P 1 + a — P 2 ' 

at oxi ax 2 opi ap 2 
d± 2 . , dx 2 . , <9x 2 . <9x 2 



x 2 =3-x 2 = - — Xl + - — x 2 + - — pi + - — p 2 . 
at oxi ax 2 opi ap 2 



(5.4a) 
(5.4b) 



Inserting appropriate PN contributions to x a and p a given in Eqs. (|4.2I) and (|4.3(l . we obtain 1PN , 2.5PN and 
3.5PN corrections to the familiar Newtonian expression for x a and it reads — apart from the 2PN and 3PN terms — 
El: 
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We now go on to present x a in terms of (x a ,v a ). As Eqs. I|5.5(l involve p a at 1PN, 2.5PN and 3.5PN orders, we use 
only Newtonian, 1PN, 2.5PN contributions to p a , given in Eqs. (|5.3[1 [rather than Eqs. (|3.8[1 truncated at the 1PN 
approximation]. After some heavy algebra, we get xi = /(xi, x 2 , vi, v 2 ) and X2 = /(xi, x 2 , vi, v 2 ) (24|: 
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In the next subsection we show the second way to com- 
pute the equations above. 



C. Second method to calculate the 
Euler-Lagrangian equations of motions 

The second route is less straightforward, but leads to 
the confirmation of the preceding result. We start from 
the Hamiltonian equations H2.3(l for p a : 



Pa 



dH 



(5.7) 



where we replace p a by the derivative of the right-hand 
side of Eq. (|5.3|) . The formal time derivatives of xi and 
x 2 , denoted by ai and a 2 respectively, are kept unevalu- 
ated. The right-hand side of Eq. H5.7|l involve the reaction 
contributions Eqs. 1)4.3(1 . as well as the Newtonian and 
post-Newtonian parts l|3.7|) . We systematically eliminate 
the particle momenta in favor of the particle coordinate 
velocities by means of Eqs. I|5.3fl (without 1/c 7 ). To cal- 
culate xi, we consider the coefficient of ai in pi, and 
invert it up to the order 1/c 7 . The equation of motion 
(|5.7|) in its present form is multiplied by the latter re- 
sult. Next, we gather all terms but ai in the right-hand 
side. Simultaneously, the acceleration X2 is isolated in 
the same way. We arrive at the coupled system: 

ai =/(xi,x 2 ,xi = vi,x 2 = v 2 ,xi = ai,x 2 = a 2 ) , 

(5.8a) 

a 2 =/(xi,x 2 ,xi = vi,x 2 = v 2 ,xi = ai,x 2 = a 2 ) . 

(5.8b) 



The various contributions are sorted in powers of 1/c, 
and labeled with two different indices. The second in- 
dex refers to the 1/c exponent, whereas the first index 
identifies the body; N stands for the Newtonian part. 



ai =Nio 
a 2 =N 20 



— fi 2 (ai,a 2 ) + — fi 5 (ai,a 2 ) + — fi 7 (ai,a 2 ) , 

c z c° c' 

(5.9a) 

-5-f 22 (ai, a 2 ) + — f 25 (ai,a 2 ) + — f 27 (ai,a 2 ) . 
c z c° c' 

(5.9b) 

The decoupling is achieved by iterative substitution of 
Eq. H5.9afl into Eq. (|5.9b(l and vice versa. We get 

ai =Ni + -^gi 2 (ai) + ig 15 (a x ) + ig 17 (ai) , (5.10a) 

a 2 =N 20 + ^g2 2 (a 2 ) + 4"g 2 5(a 2 ) + ^g2 7 (a 2 ) • 

(5.10b) 



Next, we compute ai and a2 iteratively, and insert the 
result into the right-hand sides of Eqs. (|5.10|l . The accel- 
eration has finally the following structure 



ai =/(xi,x 2 ,xi = vi,x 2 = v 2 ) , 
a 2 =/(xi,x 2 ,Xi = vi,x 2 = v 2 ) . 



(5.11a) 
(5.11b) 



We recover Eqs. (jj^ 

In order to confirm the correctness of expressions above 
for xi and x 2 , we shall check it yields the correct energy 
loss (I3.18f) . For this goal, we begin by deriving the reac- 
tion forces with the help of Eqs. (|5.6|l 



pre 



mi xi 



(5.12a) 
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F 2 cac =m 2 x 2 . (5.12b) 

After multiplying these relations by xi = vi or x 2 = v 2 
respectively, we sum to obtain 

mi (xi ■ vi) + m 2 (x 2 ■ v 2 ) = /i(xi,x 2 , vi, v 2 ) (5.13) 

This is not yet the instantaneous energy loss due par- 
ticularly to the presence of the Newtonian and post- 
Newtonian terms coming from Eqs. (|5.6I) . We then must 
search for the function g = g(xi, x 2 , vi, v 2 ) to be added 
to h as defined in Eq. Ij5.13|) in order to isolate the radia- 
tive 2.5PN and 3.5PN parts, which does not only amount 
to removing the Newtonian and post-Newtonian contri- 
butions. We shall then be able to find the average power 
by substituting to v 2 its value (|3.16|) . 

The correct energy loss is given by the time derivative 
off = i/(xi,x 2 ,pi(x a ,v a ),p 2 (x a ,v a )) truncated at the 
1PN level. The particle momenta of the 1PN Hamil- 
tonian Ij4.4|l are expressed by means of the particle co- 
ordinate velocities l|5.3|) up to the order 1/c 5 (only). 
When performing the differentiation of the resulting en- 
ergy with respect to time, xi = ai and x 2 = a 2 — com- 
ing from the derivatives of xi = vi and x 2 = v 2 — are 
replaced by expressions l|5.6l) except in the Newtonian 
term mi (ai • Vi) + m 2 (a 2 • v 2 ) where they are left un- 
evaluated. The function g — g(xi, x 2 , vi, v 2 ) is just the 
difference d£/dt — mi (ai • Vi) — m 2 (a 2 ■ v 2 ). 

By construction, h + g contains no Newtonian or post- 
Newtonian terms. It actually represents the instanta- 
neous luminosity £<g 5PN = —d£/dt, we intended to 
check, given by the negative instantaneous energy loss. 5 
The result of the averaging procedure in the case of quasi- 
elliptic orbital motion is in agreement with our prec eding 
results and with Eqs. (4.20) and (4.21) of Ref. |2^. 

VI. COMPARISM WITH IYER- WILL 
FORMALISM 

In this section, we shall show that the dissipative parts 
of Eqs. H5.6f) are compatible with the generic reaction 

5 The full result for (C^l 5pN ) is presented in Ref. |2j|. 



force of Iyer and Will depending on (2 + 6) gauge 
parameters, as the ones given in the Eqs. (4.1a)-(4.1d) 
in Ref. Q. The first two parameters at the 2.5PN order 
and the six coefficients at the 3.5PN order correspond to 
the gauge freedom. 

As a matter of fact, it is sufficient to prove that for 
certain values of the unknowns, the latter reaction force 
coincides with ours. This requires to compute the relative 
acceleration 

a = &i2 := ai - a 2 = xi - x 2 (6.1) 
from Eq. (|5.5|l for xi = /(xi, x 2 , pi, p 2 ) and x 2 = 

/(xi, x 2 , pi, p 2 ) respectively and to go to the center-of- 
mass frame 

Pi+P2=0. (6.2) 

As the accelerations we start with are function of pi and 
p 2 rather than vi and v 2 , we must simply set pi to p 
and p 2 to p in a = a(xi, x 2 , pi, p 2 ) [see Eq. (|6.1|) ]. 

The next step consists in establishing the link between 
the particle momenta p and the relative velocity v = 
V12 := vi — v 2 in order to eliminate the particle momenta 
from a = a(xi,x 2 ,p). We thus restore the Newtonian 
and post-Newtonian part in Eqs. I|4.2|) and calculate the 
relative velocity 

v = V12 := vi — V2 = xi — x 2 . (6.3) 

We make the substitutions Pi — > p and p 2 — > p in 
v = v(xi, x 2 , pi, p 2 ). The relation v = v(xi,x 2 ,p) can 
be inverted with the aid of the iterative method described 
above to get the momentum as a function of the relative 
velocity p = p(xi,x 2 ,v). Inserting p = p(xi,x 2 ,v) in 
a = a(xi,x 2 ,p) leads to a = a(xi,x 2 ,v). After intro- 
ducing the total mass M as well as the dimensionless 
mass parameter v through Eq s. I|3.10|l and [i — vM, the 
relative acceleration reads 1241: 
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with the notation 



r = (n ■ v) = (1112 • v i2 ) 



(6.5) 



This result corresponds to the generic equation of motion (1.4) and Eqs. (2) in Ref. 0] by Iyer and Will (in the 
following equations we take G = 1 = c): 



dt 2 



M M 

— r-n+ -- 



8 MM 
5"r*T 



1PN 



^2PN + ^3Pn) + rv (i?iPN + -E>2PN + £? 3 pn) 



(^2.5PN + ^3.5Pn) — V (i?2.5PN + -B3.5PN) 
I 



(6.6) 



with [sec Eqs. (2.2a), (2.2b), (2.5) and (2.12) in Ref. 0] 
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Note that we are not interested in the 2PN and 3PN 
terms in the present work. 

Iyer and Will establish in Rcfs. 0, that the energy 
and angular momentum balances hold if and only if the 
coefficients cii, &i, c% and di satisfies the following rela- 
tions: 



23 

a 1 = 3 + 3/3, 0,2 = — + 2a - 3/3, 03 = — 5/3, 



b 1 = 2 + a, b 2 = 2 - a , 63 
for the 2.5PN coefficients [see Eqs. (2.11) of Ref. @], and 
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28 2 



-2*2-5*4-6*5, 
-7*4, 
1 , * 

- ^"(1533 + 498 v) - a(14 + 9 1 

- 2 * 3 - 3 * 5 , 

d x = - 3(1 - 3 v) - ^a(l - 3 i/) - *i 



(6.9d) 
(6.9e) 

-4 1/) 

(6.9f) 

(6.9g) 



d 2 =- ^t(139 + 768 v) - -a(5 + 17 v) + *i - * 3 , 
84 2 



C?4 = 
<k 



S (369 - 624 1/) + I (3 a + 2 /?) (1 - 3 v) 
+ 3 *i — 3 *6 , 

^(295 - 335 v) + ^a(38 - 11 1/) - 3/3(1 - 3 z/) 



(6.9h) 
(6.9i) 



+ 2*4+4*3 + 3*6, (6.9j) 

=4(19 -72 1/) -5/3(1 -3 1/) + 5*6, (6.9k) 
zo 

= -^-(634- 66i/) + a(7 + 3i/) + *3 (6.91) 



for the 3.5PN coefficients [see Eqs. (2.18) of Ref. @]. 6 

The two parameters a, (3 at the 2.5PN order and 
the six coefficients Si at the 3.5PN order correspond to 
the gauge freedom and have no physical meaning. As 
mentioned in Ref. at the 2.5PN order, the values 



6 In Ref. |sj the notation £5 was used in place of i 



18 



a = — 1, 8 = characterize the gauge of Damour and 
Deruelle |2g, while a — 4, = 5 correspond to the so- 
called Burke- Thorne gauge [see for example §36.11 of [2^ ] 
also used by Blanchet 0. 

It is then a non-trivial check to verify that the 18 co- 
efficients a,i, bi, Ci and di parameterizing our 2.5PN and 
3.5PN terms [Eqs. (|6~7cjl - l|6~7f)l specialized for Eq. ijOjl ] 
yield a unique, self-consistent solution for the 8 gauge 
coefficients. They are p4|: 



for the 2.5PN order: 



a =3 



for the 3.5PN order: 



r 41 677 

6, = 1 v , 

1 84 42 ' 

. 583 157 
961 16 



= 3; 



61 151 

do = 1 

2 14 14 

r 115 15 

Oa = V , 

28 7 

51 23 

<56 = -14-14' / 



(6.10a) 

(6.10b) 
(6.10c) 
(6.10d) 



VII. DISCUSSION 

In this paper, we have computed PN reactive contri- 
butions to the ADM Hamiltonian for a pair of compact 
objects moving in general orbits, using existing results 
available for the post-Newtonian iV-body Hamiltonian, 
computed by Jaranowski and Schafer pj. The two-body 
Hamiltonian is employed to compute the first PN cor- 
rections to the instantaneous gravitational energy loss 
associated with inspiraling compact binaries. We also 
derive PN corrections to reactive equations of motion 
via Hamiltonian and Euler-Lagrangian methods. The ex- 
pressions for PN reactive accelerations, computed in the 
ADM gauge, are consistent with those of Iyer and Will. 
Employing the various reactive equations of motion we 
also computed instantaneous gravitational wave luminos- 
ity for the compact binary in general orbits. These lumi- 
nosities are in total agreement with their orbital averaged 
counterparts in the case of quasi-elliptical orbital motion. 
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APPENDIX A: HADAMARD REGULARIZATION 

Let / be a smooth real-valued function defined on 
M 3 deprived of a point Xo € K 3 , where it may be sin- 
gular. We consider the family of auxiliary functions 
/n(e) := / (xo + en), labeled by unit vectors n. We ex- 
pand / n into Laurent series around e = 0: 



/n(e) = a m (n)e r ' 

m=—N 



(Al) 



The coefficients a m of this expansion depend on the unit 
vector n. We define the regularized value of the function 
/ at Xo as the coefficient of e° in the expansion (|A1|) 
averaged over all directions: 



/rcg(x ) := <j> a (n) dQ , 



(A2) 



where d£l is the elementary solid angle viewed from the 
point xo. This procedure is called Hadamard regulariza- 
tion. To calculate the integrals in Eq. I)A2|) we used the 
following equations 



1 

-in 
1 

47T 



. n 11 dfl = , 



i 12 " dn 



for odd I . 



(A3a) 



1 



(2^+1)!!^ 



fi{ili2 _ _ _ ( J42p-li2p} 



(2p-l)!! terms 

(A3b) 



with (2 p+l) !! = 1 • 3 • 5 • . . . • (2p + 1). The curly braces 
m Eq. llA3bl mean the symmetric combination of the 
indices: 



§{ij § ki} = § ij § kl + gikgji + s u S jk 



(A4) 



3!! 



During the identification process of the primed and 
unprimed variables, xi< — * xi and X2' — > X2, in the cal- 
culation of the derivation of the Hamiltonians given by 
the Eqs. H2.6|) and (|2.8|) . some singularities appear. The 
purpose of the regularization is the elimination of these 
singularities. It is a matter of choice whether we apply 
first the regularization for index 1' and then for index 2', 
or the other way around. 

As an example, we shall apply Hadamard regulariza- 
tion procedure on the term 

(nn/ • n 22 ') (nxi/ ■ pi) (n 22 ' • Pg) 

r U> r l>2> 

= nii/n 22 ,n^ 1/ p{n 22 /P2 r ii' r i'2' • (A5) 
We intend to identify index 1 with 1' as a first step. We 



start by calculating the limit for g\ = 
replacement: 



Qi 



after the 



(A6) 
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in Eq. (|A5|) . To this purpose, we expand 



+00 



U' 



E = E ^f-tvrk'dLi® ■ (A7) 



We now perform the first angular integration over the 
terms of Eq. I|A5|I depending on the variable labeled by 
V . Here is rn = —3 and / = 2. 



1 

4-7T 



n ii' n ii' r ii' r i'2' dfiii> 



^ (—iy 1 r 

= E —jT-r l {i?d {12 , )L r- 2 3 ,— &rt lv n\ v n L lv dfi u , 
1=0 ' J 



1 

30 



By using 

d ljr m = [ m (m - 2)nV + m<S y ] r m ~ 2 . 



Ar m = m(m + l)r" 



we get 



n ll' n ll' r ll' r l'2' <^11' 



47T 



12' "12' 



12' 



(A8) 



(A9a) 
(A9b) 



(A10) 



In further computations we replace the terms on the 
right-hand side of Eq. IjAlOjl with 



' 12' 

n\ 2 , 



1 12 " 

111 
ri2> 



' 22' 



"12 + Q-2 ■ 



(Alia) 



'12 



^22' 

ri2' 



'22' 



^12 

ri2> 



'12 



Q2 

ri2> 



'22' 



(Allb) 



The second step will consist in identifying index 2 with 
2'. We calculate the limit for q\ — r 22 , — » 0. The second 
angular integration over the finite part terms (i.e. the 
terms that are of order zero in power of q\) depending 
on the variable labeled by 2' can be performed immedi- 
ately because there are no singularities in Eqs. i|All[l for 



•n\ 2 n{ 2 P{P2 



12 



The result is 



— ^n 22 ,n\ 2 , dn 7 



■n\ 2 n{ 2 P 3 iP2 



12 



-6 ik 
3 



(A12) 



(nn/ ■ n 22 ') (nir • Pi) (n 2 2' • P2) 



11" 1'2' 



(ni2 • pi) (ni2 • P2) 



°' 12 



(A13) 
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